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Abstract 

In the field-driven, thermostatted Lorentz gas the collision frequency in- 
creases with the magnitude of the applied field due to long-time correlations. 
We study this effect with computer simulations and confirm the presence of 
non-analytic terms in the field dependence of the collision frequency as pre- 
dicted by kinetic theory. 
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In a recent paper JIJ, Panja, Dorfman and van Beijeren used kinetic theory to derive 
analytic expressions for the Lyapunov exponents of the random, two-dimensional field-driven 
Lorentz gas at moderately high densities. Their approach is based on the BBGKY hierarchy 
equations and takes into account correlated collision sequences. Such "ring collisions" lead 
to long-time tails in the Green-Kubo time correlations functions affecting the transport 
properties of the system ||. Panja et al. studied the effect of ring collisions on the Lyapunov 
exponents and found that long-time correlations cause a logarithmic dependence of the 
Lyapunov exponents on the applied field strength. These non-analytic terms can be traced 
back to logarithmic terms in the field dependence of the collision frequency. Unfortunately, 
the non-analytic contributions to the Lyapunov exponents are too small to be detected with 
current computer simulation techniques jj. In contrast, the predicted effect on the collision 
frequency is large enough to be observed numerically. The purpose of this note is to verify 
the presence of logarithmic terms in the collision frequency with computer simulations. 

The random, two-dimensional field-driven Lorentz gas, shown in Fig. Jj, consists of a 
point particle of mass m and charge q moving under the influence of an external homogeneous 
field E in a two-dimensional array of circular non-overlapping scatterers with radius a. The 
scatterers are fixed at random positions in the plane and have a density of n = N/A, where TV 
is the total number of scatterers in the area A. In our numerical simulations we use periodic 
boundary conditions as indicated in Fig. [I]. When the point particle leaves the simulation 
box through one specific boundary, it re-enters it through the opposite boundary. Between 
collisions with the scatterers the point particle moves smoothly according to 
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r = — , p = gE-ap, a = q — — , (1) 
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where r = {r x ,r y } and p = {p x ,Py} are the position and the momentum of the moving 
particle, respectively. The second term in the momentum space part of the equations of 
motion is a Gaussian thermostat designed to remove the dissipated energy and keep the 
speed v = p/m of the moving particle constant [§]]. When the point particle collides with a 
scatterer it is reflected elastically, i.e., 
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v + = v_ - 2(v_ ■ u), (2) 

where u is the unit vector in the direction from the center of the scatterer to the collision 
point and v_ and v + are the pre- and post-collisional velocities of the moving particles 
respectively. 

In equilibrium, i.e. for E = 0, and for long times, the phase distribution is uniform on 
the energy shell. In this case the collision frequency of the point particle with the scatterers 
can be obtained by determining the fraction of phase space available for collision in an 
infinitesimal time interval ||: 
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In the presence of an external perturbation, however, the phase space distribution of non- 
equilibrium steady state collapses onto a multifractal strange attractor with information 
dimension strictly less than the phase space dimension. As a consequence, no analytical 
expression for the phase space distribution is available and the calculation of phase space 
integrals becomes cumbersome. Nevertheless, Panja, Dorfman and van Beijeren succeeded 
in determining the effect of the external field on the collision frequency by adopting a kinetic 
theory approach 

Qualitatively, it is evident that the external field increases the collision frequency. In the 
presence of the field, the velocity of the particle tends to be aligned in field direction before 
a collision. Immediately after the collision the velocity of the particle will therefore point 
against the field. Due to the action of the field the particle will turn around and possibly 
hit the same scatterer again after a time shorter than the average equilibrium collision time. 
The typical time scale for such a reorientation is of the order of mv / qE || . Such correlated 
collisions enhance the collision frequency. 

This argument was made quantitative by Panja, Dorfman, and van Beijeren [jl]. As a 
consequence of correlated collisions, i.e. multiple collisions of the particle with the same 
scatterer, separated by sequences of intermediate collisions, a field dependent contribution 
appears in the collision frequency: 
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where e = q\E\/(mv). The first term on the right hand side of the above equation is 
the collision frequency in equilibrium. The second, non-analytic field dependent term is 
responsible for the non-analytic dependence of the Lyapunov exponents on the field strength. 

To verify the presence of the logarithmic term in equation (Q) we have performed exten- 
sive simulations of the driven random Lorentz gas at various densities and field strengths. 
The collision frequency is obtained in a straightforward way by following the time evolution 
of the system for a long time and counting the number of collisions. For this purpose we 
use an analytical solution of equations of motion ([!]) and determine the collision point and 
the collision time of the moving particle with the scatterers numerically |fj]. Typically, we 
study systems with iV = 10 5 scatterers in a square simulation box with periodic boundary 
conditions. Only at the highest density (n = 0.15a -2 ) we use N = 5000 scatterers. By 
averaging our results over various scatterer configurations we integrate over the quenched 
disorder of the randomly placed scatterers. Computing time is saved by dividing the simu- 
lation box into cells such that at each time only a few scatterers need to be considered as 
possible collision partners. For each density and field strength we typically carry out a total 
of more than 10 9 collisions obtaining collision frequencies accurate to better than 0.01%. 
This high accuracy is needed to detect the small changes in collision frequency in the weak 
field regime. 

For analysis it is convenient to rewrite equation (f|) as 

8v a 2 a ea 

— = m(2n<r) In — (5) 
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where 5v is the deviation of the collision frequency from its equilibrium value. Accordingly, 
bvje 2 should behave linearly when plotted as a function of ln(ea/v). Figure ^| shows 8vje 2 as 
a function of ln(ea/v ) for the densities n = 0.001a" 2 , 0.002a" 2 , 0.005a" 2 , 0.01a" 2 , 0.02a" 2 , 
0.04a" 2 , 0.08a~ 2 , and 0.15a~ 2 . In Fig. ||] densities increase from left to right. Each data 
point is obtained as an average over 50 to 500 runs with different scatterer configurations. 



The error bars are estimated from the variation of the collision frequency in these sets of 
runs. The solid lines connect data points corresponding to the same density and the dotted 
lines are straight lines with slope —a/(27rv) fitted to the data in the low field regime. 

For all densities but the highest one Su/e 2 is a linear function of ln£ in the low field 
strength range with slope —a/(2nv) in agreement with equation @. This confirms the 
existence of the non-analytic field dependent term predicted by kinetic theory. 

From Fig. |2| we can also infer the range of validity of the theory in the n-e plane. 
Clearly, equation (||) breaks down for densities larger than about n = 0.08a -2 . For densities 
below n = 0.08a -2 equation (|5|) holds for field strengths below a certain critical value e c . 
The intercept of the fitted lines with the x-axis can be used as a measure of this critical 
value. e c is an approximately linear function of the density growing from e c ~ 0.006t>/a at 
n = 0.001a" 2 to e c ~ 0.5v/a at n = 0.08a" 2 . 

However, the agreement between numerical results and analytical theory is not perfect 
even for low densities and weak external fields. While the lines corresponding to equation 
(H) have the correct slope, their intercept with the y-axis is lower by about 0.2a/ v than the 
value obtained numerically. This discrepancy is essentially constant in the density range we 
have studied. It must be due to the approximation of the probability of return to a given 
scatterer by a solution of the diffusion equation, as described in the Appendix of Ref. [|TJ. 
In order to check this in detail one would have to explicitly solve the Lorentz-Boltzmann 
equation in presence of the field and thermostat. 

In summary, we have unequivocally detected the appearance of non-analytic terms in 
the collision frequency and hence, indirectly, in the Lyapunov exponents of the random 
field driven Lorentz gas. Our numerical results agree quantitatively with the kinetic theory 
predictions of Panja, Dorfman, and van Beijeren. 
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FIGURES 

FIG. 1. Geometry of the random, field-driven Lorentz gas. 

FIG. 2. bvje 2 as a function of ln(ea/v) for the densities n = 0.001a" 2 , 0.002a' 2 , 0.005a" 2 , 
0.01a" 2 , 0.02a" 2 , 0.04a" 2 , 0.08a" 2 , and 0.15a" 2 (from left to right). The dotted lines are straight 
lines with slope —a/{2'Kv) fitted to the data. The linear behavior of 5v/e 2 in the weak field regime 
indicates the presence of the logarithmic terms predicted by kinetic theory. 
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Fig. 1 Dellago, v. Beijeren, Panja, Dorfman 
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Fig. 2 Dellago, v. Beijeren, Panja, Dorfman 



